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Abstract 

The results of a new search for positronium decays into invisible final states are 
reported. Convincing detection of this decay mode would be a strong evidence 
for new physics beyond the Standard Model (SM): for example the existence of 
extra-dimensions, of milli-charged particles, of new light gauge bosons or of mirror 
particles. Mirror matter could be a relevant dark matter candidate. 

In this paper the setup and the results of a new experiment are presented. In a 
collected sample of about (6.31 ± 0.28) x 10 6 orthopositronium decays, no evidence 
for invisible decays in an energy window [0,80] keV was found and an upper limit 
on the branching ratio of orthopositronium o — Ps — > invisible could be set: 

Br(p -Ps^ invisible) < 4.2 x 10~ 7 (90% C.L.) 



Our results provide a limit on the photon mirror-photon mixing strength e < 
1.55 x 10~ 7 (90% C.L.) and rule out particles lighter than the electron mass with 
a fraction Q x < 3.4 x 10~ 5 of the electron charge. Furthermore, upper limits on 
the branching ratios for the decay of parapositronium Br(p — Ps — > invisible) < 
4.3 x 10 -7 (90% C.L.) and the direct annihilation Br(e + e~ — ► invisible) < 2.1 x 
10~ 8 (90% C.L.) could be set. 
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1 Introduction 



Although direct manifestations of new physics are searched for at the high 
energy frontier, new phenomena can also be looked for at low energies via 
precision measurements. The new effects might be observed in rare decays of 
the positronium (see e.g. [1]). 

Positronium (Ps), the positron-electron bound state, is the lightest known 
atom, which at the current level of experimental and theoretical precision is 
bound and self-annihilates through the electromagnetic interaction [2]. This 
feature has made positronium an ideal system for testing the accuracy of 
Quantum Electrodynamics (QED) calculations for bound states, in particu- 
lar for the triplet (l 3 ^) state of Ps (orthopositronium, o-Ps) [3]. Due to the 
odd-parity under C-transformation o-Ps decays predominantly into three pho- 
tons with a lifetime in vacuum of r Q _p s = 142.05 ns [2]-[5]. The singlet (l 1 ^) 
state (parapositronium, p-Ps) decays predominantly into two photons with a 
lifetime in vacuum of r p _ Ps = 125 ps [6,7]. The longer lifetime of o-Ps (due to 
the phase-space and additional a suppression factors) gives an enhancement 
factor ~ 10 3 in the sensitivity to an admixture of potential new interactions 
not accommodated in the Standard Model (SM) [9] . 

This paper focuses on a new search for o — Ps — > invisible decays. By invisible 
we mean photonless decays, i.e. decays which are not accompanied by energy 
deposition in a hermetic calorimeter. In the SM the decay into a neutrino- 
antineutrino pair has a branching ratio of 6.6 x 10~ 18 [8]. Evidence for invisible 
decays in the region ~ 10 -7 would, therefore, unambiguously signal the pres- 
ence of new physics. New models that are relevant to the o — Ps — >• invisible 
decay mode predict the existence either of i) extra-dimensions [10,11], or ii) 
fractionally charged particles [12,13], or iii) a new light vector gauge boson [14], 
or iv) mirror particles, which could be candidates for dark matter [15]-[20]. 

The first experiment to search for invisible decay channels of o-Ps was per- 
formed by Atoyan et al. [21]. Their result on Br(o — Ps — > invisible) < 
5.3 x 10~ 4 (90% C.L.) excluded this channel as a possible explanation of the 
o-Ps lifetime anomaly (for a recent review see e.g. [22]). This search was re- 
peated by Mitsui et al. who found a branching ratio Brio — Ps — > invisible) 
< 2.8 x 10~ 6 (90% C.L.) [23]. Furthermore, they could place a limit on the 
existence of milli-charged particles and on the photon mirror-photon mixing. 
This result was corrected in [24] by taking into account the suppression factor 
for the mixing due to the presence of matter. 

The rest of the paper is organized as follows. The details of the experimental 
setup are reported in Section 2. The expected background is presented in 
Section 3. In Sections 4 and 5 the data analysis and the results are described. 
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The interpretation and the conclusion are reported in Sections 6 and 7. 



2 Experimental technique and setup 

The experimental signature of o — Ps — >■ invisible decay is the apparent 
disappearence of the energy 2m e expected in ordinary decays in a hermetic 
calorimeter surrounding the o-Ps formation target. The readout trigger for the 
calorimeter is performed by tagging the stopping of a positron in the target 
with high efficiency 

For the design of the experimental setup aiming at a sensitivity Br(o — Ps — > 
invisible) ~ 10~ 8 , the following criteria were considered: 

(a) The probability not to detect all direct e + e~ annihilation photons was 
suppressed to < 10~ 9 using a thick hermetic crystal calorimeter, minimiz- 
ing the dead material and with a good energy resolution. The probability 
to loose all photons in 37 decays is consequently even smaller. 

(b) The region around the target was designed with as little dead material 
as possible in order to reduce photon losses. 

(c) By the appropriate choice of a porous target material and its dimensions 
a high fraction of o-Ps was produced resulting in a suppression of the 
background from the the direct e + e~ annihilation and p-Ps decays and 
in high statistics. 

(d) The trigger rate and the DAQ speed were maximized for statistics. 

(e) An efficient positron tagging system was designed to provide a clean 
trigger for positronium formation. A method was developed to suppress 
the background from the electrons emitted in the EC process (shake-off 
electrons) . 

(f) An efficient identification of the 1.27 MeV photon emitted by the 22 Na 
radioactive source was achieved with a method to veto the charged par- 
ticles entering the trigger counter, thus, the backgrounds related to them 
could be reduced. 

The schematic illustration of the detector setup is shown in Figure 1 (the 
detailed description of the experimental technique and setup can be found in 
Refs. [25,26]). 

Positrons are produced from a 22 Na source with an activity of ~ 30 kBq. The 
22 Na has a half life of 2.6 years and has a Q-value for the nuclear transition 
to 21 Ne of Q = 2.842 MeV. This is the maximum energy available for the 
particles involved in one of the three possible decay modes of 22 Na: 

(a) Decay mode A (Br ~ 90.6%): the (3 + decay with end-point energy 0.546 MeV. 
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Fig. 1. Schematic illustration of the experimental setup: a) front view, b) top view. 

The positron is always followed by the prompt emission of a 1.27 MeV 
photon (r ~ 3.7 ps) from the 21 Ne* de-excitation to the ground state. 

(b) Decay mode B (Br ~ 9.44%): the Electron Capture process (EC), where 
an orbital electron is captured by the nucleus, and only a 1.27 MeV 
photon and a neutrino are emitted from the source. In some rare cases 
(with a probability ~6x 1CT 3 ), an orbital electron is ejected, due to the 
sudden change in the nucleus charge (shake-off) [28]. 

(c) Decay mode C (Br ~ 0.056%): there is no photon emission because the 
transition goes directly to the ground state. The end point energy of the 
positron is 1.83 MeV. 

Therefore, in most o-Ps (p-Ps) decays we expect 3(2) photons with a summed 
energy equal to 2m e and one photon with an energy of 1.27 MeV (see Figure 2). 
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Fig. 2. Schematic illustration of the positron tagging system and the o-Ps formation 
target of the setup. 
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Photons from the direct e + e~ annihilation in flight or from the positronium 
decays were detected in a hermetic, segmented BGO calorimeter (the ECAL). 
Two endcap counters called TBGO and FBGO (see Figures 1-2) surrounded 
the target on each side. At the analysis level the 1.27 MeV photon ("the 
triggering photon") was required to be identified in the TBGO counter. The 
calorimeter was instrumented with charge and time readout. 

The activity of the source was chosen to maximize the trigger rate versus the 
inefficiency of signal detection (mostly due to pileup events). The source was 
prepared by evaporating drops of a 22 Na solution directly on a 100 fim thick 
and 2x8 mm 2 wide plastic scintillator (see Figure 3) fabricated by squeezing a 
500 yum diameter scintillating fiber (Bicron BF-12). In this way, no dead mate- 
rial was introduced for a source holder. The S-shape of the fiber (see Figure 1) 
was selected to avoid background from back-to-back 511 keV annihilation pho- 
tons. The scintillating fiber was read out at both ends by two photomultipliers 



Fig. 3. Schematic view of the scintillating fiber with the Na source on the squeezed 
part in the center. 

(Philips XP2020) located outside the detector (see Figure 1). The coincidence 
of the two PMT signals was used to tag the passage of a positron through the 
fiber and acted as a start signal for the data readout system. The use of two 
PMTs in coincidence, instead of a single one [23], lowered the ratio between 
fake and real positron triggers to < 1.9 x 10~ 10 . 

Opposite to the source, a 4x8x8 mm 3 Si0 2 aerogel piece (type SP30, purchased 
from Matsushita Electric Works) was placed in contact with the squeezed 
fiber (see Figure 2 and 4). Positrons stopping in the aerogel target may form 
positronium (the formation probability is 45% [29]) which can migrate into the 
aerogel pores. The collisions with the walls of the pores did not appreciably 
quench the o-Ps: when the aerogel was flushed with nitrogen, a fit to the 
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distribution of the time difference between the start from the fiber and the 
stop from the calorimeter yielded r G _p s = 129.1 ±1.8 ns, which is very close to 
the lifetime in vacuum r _ Ps ~ 142.05 ± 0.02 ns [30]. 

The ECAL was composed of 100 BGO crystals that surrounded the target 
region providing a nearly isotropic sphere of radius 200-220 mm (see Figure 1). 
Each crystal had a hexagonal cross-section of 61 mm diameter and a length 
of 200 mm. The crystals of the inner most ring were wrapped in a 2 /iin thick 
aluminized Mylar foil to minimize the photon energy absorption. The other 
crystals were wrapped with Teflon foils of 750 fxm thickness. The crystals in 
the barrel and the FBGO were readout with ETL 9954 photomultiplier tubes. 

The energy deposited in the fiber is an important parameter in order to reject 
the background from the electrons emitted in the EC process (shake-off elec- 
trons, see Section 3). The mean number of photoelectrons detected in each 
XP2020 for a positron crossing the fiber was measured to be about 1.2, thus, 
a cut on the energy deposited in the fiber using these signals would not be 
meaningful. For this reason the FBGO was also used to measure the energy 
deposited in the fiber by the positron: the light emitted by the fiber could 
traverse the transparent aerogel and enter the FBGO through an aperture in 
the wrapping on its front face. This light was then guided by the FBGO to 
the PMT attached on the back face of FBGO (as illustrated in Figure 4). 
This method provided a mean number of photoelectrons equal to 13±1 for a 
positron traversing the fiber [26]. 

The TBGO endcap was used in the off-line analysis to identify the 1.27 MeV 
photon. The energy resolution of the TBGO was an essential parameter to 
reduce backgrounds related to the misidentification of the 1.27 MeV photon 
(see also Section 3). To provide a better energy resolution, this crystal was 
coupled to an ETL 9964 PMT with a more uniform light collection and a larger 
quantum efficiency than the ETL 9954. For the same reason this crystal was 
of a better quality than the others and efforts were dedicated to maximize 
the light collection by keeping the amount of dead material introduced by the 
crystal wrapping as small as possible. The best results were achieved with the 
crystal wrapped in the 3M radiant mirror (64 /xm thickness): the resolution 
at 662 keV was measured to be about 15% (FWHM). To select the triggering 
photon, an energy window [1275 ± 67] keV was used in the analysis. In addi- 
tion, to veto charged particles (positrons and electrons) entering the TBGO, 
a 1 mm thick plastic scintillator (Bicron BC-400) was optically coupled to 
the TBGO front face (for more details see [27]), i.e. the same PMT was used 
to detect the light signals from the plastic scintillator and the TBGO crystal 
(see Section 3 for a discussion of backgrounds associated to charged particles 
entering the TBGO). The signals from the plastic scintillator and the BGO 
could be distinguished because of the different decay time (r = 2.7 ns for the 
plastic scintillator and r ~ 300 ns for the BGO). For this purpose, split signals 



6 




Fig. 4. Schematic illustration of the method to readout the energy deposition in the 
fiber. 

from the PMT were fed into two ADCs with different integration gates (called 
short gate and long gate). 

A VME system interfaced to a PC was used for data acquisition (the DAQ 
rate was about 1800 events/s). For every trigger, five CAEN 32 channels QDC 
v792 modules recorded the charge of the crystal signals while a CAEN TDC 
v775 recorded the time information. A trigger gate length of 2.9 (is was chosen 
to keep the probability of o-Ps to decay after this time to < 10~ 9 . 

A temperature stabilized light-tight box containing the calorimeter was built 
in order to keep the temperature of the BGO crystals in the range of ±0.5 °C. 
Water, whose temperature was controlled by two thermostats, circulated through 
copper tubes welded on two copper plates inside the box. The experimental 
hall was air conditioned to keep temperature variations within ±1°C. The 
high-voltage dividers of all PMTs were placed outside the box in order to avoid 
energy dissipation close to the crystals. The BGO crystals were equipped with 
LEDs that could be pulsed periodically to monitor the response. Additionally, 
the gains of the PMTs were also monitored to check their stability. 

The detector was calibrated and monitored internally using the 511 keV an- 
nihilation photon and the 1.27 MeV photons emitted by the 22 Na source. 
Variations of the energy scale during the run period were within < 1% and 
corrected on the basis of an internal calibration procedure [26]. 
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3 Background estimation and dedicated engineering run 



In order to reach the required sensitivity, the background must be reduced 
and controlled at the level of 1CT 8 . To understand the different background 
sources and to cross-check the simulation, we perfomed an engineering run 
with a simplified version of our detector. During two months of data taking, the 
stability of the detector and its components was investigated. The comparison 
between the background expected based on Monte Carlo (MC) simulations 
and the data of the engineering run is summarized in Table 1. 

The ECAL thickness of 200 mm provides a probability of < 10~ 9 for two 
511 keV photons to escape detection (see background 1 in Table 1). For three 
photons decays, this probability is consequently even smaller. 

If one (or more) annihilation photon (e.g. backscattered from the target) over- 
laps with the 1.27 MeV in the TBGO, it can fall in the trigger energy window 
[1275 ± 67] keV because of the finite energy resolution. This introduces a 
background if the remaining annihilation photon gets absorbed in the dead 
material or escapes detection. The separation between the upper bound of the 
trigger energy window and the sum of a triggering and a 511 keV photon was, 
thanks to the good energy resolution of the TBGO, 7a of the 1786 keV peak 
(1275+511 keV). Thus, the level of this background is < 5 x 10~ 9 (background 
2 in Table 1). 

In order to suppress the other sources of background related to the misidentifi- 
cation of the 1.27 MeV photon (backgrounds 3, 4 and 5 listed in Table 1), one 
had to veto charged particles entering the TBGO. Two processes are respon- 
sible for generating such triggers. One is associated with decay mode B (EC) 
when the 1.27 MeV photon interacts in the fiber, faking a positron signal. If 
the scattered photon and the Compton electron reach the TBGO, the sum of 
the energy of the two particles can be misidentified as the triggering photon 
without any energy in the rest of the calorimeter. Another background can 
occur in decay mode A if the 1.27 MeV photon is not detected: a trigger can be 
produced by a positron that multiple scatters (MS) in the fiber and deposits 
enough energy to trigger the experiment. If the positron reaches the TBGO 
with a kinetic energy of about 200-300 keV and the two 511 keV annihilation 
photons are completely absorbed in the TBGO an energy close to 1.27 MeV 
will be reconstructed. This will appear as an invisible decay since no energy 
is expected in the rest of the detector. If the 1.27 MeV photon is not present 
(decay mode C) a trigger can similarly be produced by the 1.83 MeV positron. 
To veto these backgrounds the charged particle veto of the TBGO described 
in the previous section was used such that the backgrounds 3, 4 and 5 had a 
probability of < 10~ 8 in the final setup. 
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The EC photon may accidentally coincide with a trigger from the fiber gener- 
ated either by the PMTs noise or by some other particles emitted from possible 
unstable isotopes formed during the target activation. The level of this back- 
ground could be reduced with the selection of two XP 2020 with very low 
noise (< 30 counts/s) and the requirement of the coincidence between them. 
In addition, a radioactive source with a controlled high purity was chosen. 
From the data of the engineering run, this background was estimated to be 
< 1.9 x 10~ 10 (background 6 in Table 1). 

During the decay mode B (EC) the fiber signal can be generated by shake- 
off electrons (background 7 in Table 1). Since the probability of an electron 
ejection steeply drops with its emission kinetic energy (more than 4 orders of 
magnitude in the first 100 keV) the cut on the energy deposited in the fiber 
by the triggering particles was used to suppress this background. 

The engineering run allowed to test these backgrounds as shown in Table 1. 
The expected fraction of zero energy events is 10% smaller than what was mea- 
sured. This difference can be explained by the contribution from the shake-off 
electrons which was not included in the simulation. Indeed, in the engineer- 
ing run there was no information about the energy deposited in the fiber so 
that no cut on the energy of the particles passing through the fiber could be 
applied. 

The last column of the table lists the expectations for the final setup. The total 
background is estimated to be at the level 10~ 8 . The threshold on the energy 
deposited in the fiber was set considering the uncertainty due to the number of 
photoelectrons and the subtraction method to be sure that no electron below 
100 keV could trigger the fiber; it was chosen, based on simulations, to be at 
140 keV. 



4 Data analysis 

For the analysis we used a data sample of 1.39 x 10 10 recorded fiber triggers 
collected over a four months data taking period. For each event the following 
variables were used and cuts were applied to suppress background: 

(1) AT short is the time from the trigger start to the end marker of the dual 
timer unit that generates the short gate (measuring the light from the 
plastic scintillator coupled to the TBGO). 

(2) Xiong is the pedestal of one of the QDC channels integrated using the long 
gate as a start trigger. 

(3) ATxp is the time difference between the two XPs reading the fiber. 

(4) AT TBGO is the time difference between the XPs coincidence and the 
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Table 1 



Comparison between expected and measured background level for the different back- 
ground sources in the engineering run and expected background level for the final 
setup (see text for details). 

TBGO. 

(5) -Stbgoc is the energy deposited in the TBGO with a gate delayed by 
15 ns to measure the energy deposited in the TBGO crystal without the 
contribution of the plastic scintillator. 

(6) -Stbgo is the energy deposited in the TBGO with the full gate (long 
gate). 

(7) -Efbgo is the energy deposition in the fiber measured with the FBGO. 

The distributions of these variables for a reduced data sample of 10 6 triggers 
are shown in Figure 5. The used selection cuts are listed in Table 2. The cuts 
were selected and tuned with the help of a dedicated run with corresponding 
statistics of about 5% of the data. These variables can be grouped in three 
categories, depending on their function: (a) The first two variables check the 
stability of the electronics and the duration of the gate widths. The selection 
has been tuned experimentally looking at the obtained spectra, (b) The vari- 
ables 3) and 4) suppress accidental triggers faking positrons in the fiber. The 
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Fig. 5. The cuts applied to the variables. The numbers on the plot correspond to the 
variable defined in the text. Only the colored regions contribute. 



cuts have been chosen to minimize the accidentals and maximize the signal 
statistics, (c) The variables 5) to 7) are the cuts that reduce triggers from 
backgrounds that mimic the appearance of a positron in the formation cav- 
ity region. Furthermore, the upper limit of the energy window for the 1.27 
MeV photon is very sensitive to the background from the "absorption" of one 
511 keV 7 in the trigger energy window. Therefore, the long gate integration, 
which has a better resolution, is used to define this selection, while the short 
gate is used to reject the events with some energy deposited in the scintillator. 
The cut of ±lcr was selected in order to enhance good triggers to the required 
level. 

All the selection cuts, except 7), have been defined in terms of the sigma of 
the signal determined with a Gaussian fit to the data sample with 10 6 events. 
Table 2 summarizes the values used and the evolution of the total number of 
events passing the cuts. 
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The lower cut for the energy deposited in the fiber has been chosen to reduce 
the probability of shake-off electrons to trigger the fiber to a level < 10~ 8 as 
discussed in Section 3. The measured fraction of the o-Ps produced in the 
aerogel is reduced by 20% applying the threshold for the energy cut in the 
fiber. This was expected, since the positrons that deposit the most energy are 
the ones stopping in the fiber. 



Variable 




Selection cut 


Fraction of events 


name 


#a's 


value of la 


remaining after cuts 


1) AT s hort 
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3.71 ns 
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5) Etbgo 


±la 


74 keV 


3.1% 


6) E TBGOc 


±1(7 


67keV 


2.7% 


7) EpBGO 


140 keV < E FBGO < 400 keV 


1.1% 



Table 2 

Definition of cuts and the remaining fraction of events after the cut is applied. 



5 Results 

After imposing the above requirements a final sample of 1.41 xlO 8 events was 
obtained. For these events the energies of all the 100 BGO crystals, except the 
TBGO, were summed. Figure 6 shows the spectrum of the total energy (E tot ) 
deposited in the ECAL. The peak at 1022 keV corresponds to the positronium 
mass (M Ps ~ 2m e ). The inset shows that no event is observed in the zero 
energy region. 

To define the upper energy cut on E tot , below which an event is considered as 
photonless (invisible), a dedicated run of 10 7 triggers was performed triggering 
the experiment only with the 1.27 MeV photon and no requirement of the 
fiber [21]. The zero peak contained about 10.4±1.2% of the events passing 
the selection cuts defined in Section 4. This value was corrected by a factor 
0.8 determined using the MC simulations to take into account the different 
detection efficiencies of the 1.27 MeV photon in the case of the EC process 
and in the case of the transition with the positron. The measured value was 
consistent with the expected fraction of electron capture events (decay mode 
B). The cut E tot < 80 keV corresponding to the region containing 99% of the 
events in the EC peak at zero energy was used to define the photonless events. 
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To determine the signal inefficiency, mostly due to pileup, ~ 10 6 events have 
been collected using a random trigger formed by delaying the fiber coincidence 
by 16 /is (half of the mean time interval between two events). For the 80 
keV threshold defined above, this gave an inefficiency of (11.6 ± 0.5)% that 
was consistent with the prediction of the simulation. This inefficiency was 
measured at the beginning of the data taking and, conservatively, was not 
corrected for the reduction of the source intensity during that period. 



For the FBGO, the broadening of the pedestal due to the contribution of 
energy deposited in the fiber had to be considered. The energy of the two 511 
keV annihilation photons was localized in two opposite-located crystals of the 
barrel and then the FBGO pedestal was built. This resulted in a correction of 
1% in the efficiency of zero signal detection. 



Finally, the signal efficiency was estimated to be e ~ (87.4 ± 0.5)%. 




Fig. 6. Spectrum of the sum of the total energy in the ECAL. The inset shows the 
magnified view of the low-energy region in logarithmic scale. 



The mean fraction of o-Ps in the data sample could be evaluated from the de- 
cay time curve by fitting the observed distribution to the function A-e ( -~ t ^ T °- p ^ + 
B (B is the accidental background) starting from the time t = 100 ns when 
o-Ps was completely thermalized in the target [26]. After taking into account 
the estimated difference of efficiency for 2 and 3 gamma detection and the 
pick-off effect (measured from the lifetime spectra with the same method as 
described in Ref. [31]), the fraction of o-Ps in the data sample was 4.5 ± 0.2 % 
[26] . Since in the signal region no zero energy events were observed the upper 
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limit for the branching ratio [30] is: 



Br[o -Ps^ invisible) = 2.3/(N _ Ps ■ e) < 4.2 x KT 7 (90% C.L.) (1) 

where iV _p s = (6.31 ±0.28) x 10 6 is the number of o-Ps in the selected sample. 
Figure 6 shows the extrapolation into the region of the zero signal with an 
exponential. The integral from to 80 keV of the function obtained from the 
fit gives an evaluation of the background contribution in this region. The result 
is iVbkg = 0.34 ± 0.04 expected events, where the error was evaluated from the 
uncertainty related to the extrapolation procedure itself. 

This experiment can also be used to obtain upper limits on Br (p-Ps — > 
invisible) and on Br(e + e~ — > invisible). For this purpose the different prob- 
abilities of positrons stopping in the fiber and in the aerogel were calculated 
with the help of simulations. Several papers were reviewed in Ref. [23] conclud- 
ing that the probability to form positronium (75% of o-Ps and 25% of p-Ps) 
per positron stopping is about 0.45 in the aerogel and between 0.2 and 0.4 
in the fiber. In the aerogel a fraction of 0.05 to 0.1 of the p-Ps atoms experi- 
ence pick-off annihilation in about 2 ns before escaping out of the silica grains 
in the pores. For o-Ps this probability ranges from 0.28 to 0.45. In a plastic 
scintillator (the fiber) the corresponding pick-off probabilities are 0.99 to 1 
for o-Ps and 0.05 to 0.1 for p-Ps [23]. From the simulations, the fraction of 
positrons stopping in the fiber is 0.43 and in the aerogel is 0.25, therefore, the 
smallest possible fraction of p-Ps decays (including pick-off) is 5.5%. Thus, an 
upper bound for the invisible decay of p-Ps can be calculated: 

Br(p -Ps^ invisible) = 2.3/(N p _ Ps • e) < 4.3 x 10~ 7 (90% C.L.) (2) 

where N p _ Ps = 0.055 x 1.41 x 10 s ~ 6.14 x 10 6 was used and the positrons that 
do not stop in the fiber or in the aerogel are assumed to annihilate directly. 
The number of e + e~ decays can be calculated subtracting from the total num- 
ber of events the p-Ps and o-Ps decays, thus, one obtains N e + e ~ ~ 1.29 x 10 s 
and an upper limit for the branching ratio of: 

Br(e + e- -> invisible) = 2.3/(N e+e - ■ e) < 2.1 x 10~ 8 (90% C.L.) (3) 



6 Interpretation 

No event consistent with an invisible decay was found in the large sample of 
events. 
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Using Eq. (32) of Ref.[14], the bound for particles with a fraction Q x of the 
electron charge can be plotted as a function of their mass mx for mx < Tn e , as 
shown in Fig. 7a. Thus, the region of the charge-mass parameter space, which 
was not excluded directly by the SLAC results [32] and the previous search 
for o — Ps — > invisible [23], is covered by this experiment (see Fig. 7b). 
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Fig. 7. a) Mass-charge parameter space for the o-Ps decay into milli- charged parti- 
cles, excluded with this experiment, b ) Comparison of our results ( the dashed region 
on the plot) with other experimental (SLAC [32] and previous o — Ps — > invisible 
[23]) and astrophysical bounds (the plot was taken from [13]). 

The strength of the photon mirror-photon mixing e can be extracted from the 
limit on the Br(o — Ps — > invisible) with [24]: 



Br(o — Ps — > invisible) FsmFcoII 
2T/^ 2(1 - Br(o - Ps -»• invisible)) 



(4) 



by substituting a conservative value of r co n = 5 x 10 4 s _1 for the collision 
rate of the o-Ps against the walls of the aerogel pores [17]. Tsm is the decay 
rate of o-Ps in vacuum and / = 8.7 x 10 4 MHz is the contribution to the 
ortho-para splitting from the one-photon annihilation diagram involving o-Ps 
[16]. Using our result one can estimate the mixing strength to be e < 1.55 x 
1(T 7 (90% C.L.). This is close to the BBN limit of e < 3 x 1(T 8 [33] but does 
not cover all the region of interest suggested by the DAMA and CRESST 
results [34] and motivated by GUT predictions [18] and by string theory [35] 
(e > 10 "). 



7 Conclusion 



In this paper the results of a new search for an invisible decay of o-Ps were 
reported. Since no event was found in the energy window [0,80] keV, an upper 
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limit for the branching ratio was set: 

Br[o - Ps -> invisible) < 4.2 x KT 7 (90% C.L.) (5) 
improving the best existing bound [23] by a factor 7. 

Analyzed in the context of theoretical models the negative result provides 
an upper limit on the photon mirror-photon mixing strength e < 1.55 x 
10~ 7 (90% C.L.) and rules out particles with a fraction Q x < 3.4 x 10~ 5 (for 
rax < m e ) of the electron charge (milli-charged particles). Furthermore, an 
upper limit on the branching ratios for the process Br(p — Ps — > invisible) < 
4.3 x 10- 7 (90% C.L.) and Br(e + e~ -> invisible) < 2.1 x 10~ 8 (90% C.L.) 
could be set. 
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